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We present a comprehensive study of the optical properties of InAs/InP self-assembled quantum 
dots (QDs) using an empirical pseudopotential method and configuration interaction treatment of 
the many-particle effects. The results are compared to those of InAs/GaAs QDs. The main results 
are: (i) The alignment of emission lines of neutral exciton, charged exciton and biexciton in InAs/InP 
QDs is quite different from that in InAs/GaAs QDs. (ii) The hidden correlation in InAs/InP QDs is 
0.7 - 0.9 meV, smaller than that in InAs/GaAs QDs. (iii) The radiative lifetimes of neutral exciton, 
charged exciton and biexciton in InAs/InP QDs are about twice longer than those in InAs/GaAs 
QDs. (v) The phase diagrams of few electrons and holes in InAs/InP QDs differ greatly from those 
in InAs/GaAs QDs. The filling orders of electrons and holes are shown to obey the Hund's rule and 
Aufbau principle, and therefore the photoluminescence spectra of highly charged excitons are very 
different from those of InAs/GaAs QDs. 

PACS numbers: 73.21. La, 71.35.Pq, 78.67.Hc 



I. INTRODUCTION 

Self-assembled quantum dots (QDs) have attracted 
great interest in the past two decades due to their po- 
tential applications in the optoelectronic devices 1 - - — and 
quantum information science**— , such as high-efficiency 
and low-threshold lasers^ single photon source j2r— en- 
tangled photon emitters ; 11 ' 12 and qubits etc— These 
applications all take the advantages of the unique prop- 
erties of the QDs which have discrete energy levels and 
sharp absorption/emission finest— due to the three- 
dimensional confinement effects. 

Until very recently, most of the experimental and the- 
oretical works focus on InAs/GaAs QDs 1 -' 4 -'^— . Huge 
progress has been made both in understanding of the 
fundamental physics and towards real applications of 
these QDs. For examples, the entangled photon sources 
from biexciton cascade process have been demonstrated 
experimentally ] 11 ' 12 The strong coupling between QDs 
and cavity has been achieved— 2r— The coherently manip- 
ulation of single charge and spin in InAs/GaAs QDs have 
also been demonstrated— >2Ir— These achievements make 
QDs extraordinary promising candidates for quantum 
information applications. Theoretically, the electronic 
structures and optical properties of InAs/GaAs QDs have 
been studied intensively via both k • p methods^ - — and 
the microscopic models! 25 ' 41 ' 42 

On the other hand, InAs/InP QDs have attracted in- 
terest only very recentlyi^ril mainly because their emis- 
sion wavelengths are naturally around 1.55 /im (C-band), 
which are ideal for fiber optical telecommunication ap- 
plications. Compared to InAs/GaAs QDs, there are 
much fewer works on InAs/InP QDs, both experimen- 
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tally and theoretically. It has been shown from atomistic 
pseudopotential calculations that the electronic struc- 
tures of InAs/InP QDs differs greatly from those of 
InAs/GaAs QDs^ even though they have the same dot 
material, but different matrix material. Interestingly, it 
has been found 4 ^ that the fine structure splittings (FSS) 
in InAs/InP(001) QDs are very small, which is very suit- 
able for C-band entangled photon source. Indeed, recent 
experiments indicate that InAs/InP QDs do not have sig- 
nificant FSS ,£2 

However, the optical properties of InAs/InP QDs has 
not yet been systematically studied. In this work, we 
present a comprehensive study of the optical properties of 
InAs/InP QDs using an atomistic empirical pseudopoten- 
tial method and the configuration interaction (CI) treat- 
ment of the many-particle effects— - The methods have 
been successfully applied to InAs/GaAs QDs, and ob- 
tained many results that are very good agreement with 
experiments i 5 1 1 52 A nice review of the methods can be 
found in Ref. [53|. We find that the photolumines- 
cence (PL) spectra of InAs/InP QDs are very different 
from those of InAs/GaAs QDs in many aspects, includ- 
ing the alignment of emission lines of the neutral exciton, 
charged excitons and biexciton, the hidden correlation, 
the radiative lifetimes, and the highly charged PL spec- 
tra. We hope this work can provide helpful insight for 
the optical properties of InAs/InP QDs. 

The rest of the work is organized as follows. In Sec. HH 
we outline the methods used in our calculations. In Sec. 
IIII Al and Sec. HUB I we discuss the binding energies and 
the hidden correlation in InAs/InP QDs. The radiative 
lifetimes of exciton, biexciton, and trions are discussed in 
Sec. IIII CI The phase diagrams of electrons and holes in 
InAs/InP QDs are presented in Sec. IIII D[ whereas the 
highly charged exciton PL spectra in InAs/InP QDs are 
discussed in Sec. IIII El We summarize in Sec. IIVI 
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II. METHODS 

Figure Q] depict the geometry of the QD system used in 
our calculations. The lens-shaped InAs QD is embedded 
in the center of the 60 x 60x60 InP 8-atom unit cell on 
the top of a two-monolayer wetting layer. The single- 
particle electronic structures of the QDs are calculated 
via an empirical pseudopotential method , 25 i 41 whereas 
the many-particle effects are treated via a configuration 
interaction method, as follows. 

We first obtain the atomic positions of the QDs struc- 
ture. This is done by minimizing the strain energy of 
the system via the valance force field (VFF) method ) 54 i 55 
which has been demonstrated to be a good approxi- 
mation for semiconductors. Once we obtain the re- 
laxed positions Rj jC « of ah the atoms of type a at site 
i, we calculate the total electron- ion potential V ps (r) — 
Vso + SiS a u Q (r — Hi, a ) as a superposition of the local 
screened atomic pseudopotentials v a (r), the total (non- 
local) spin-orbit (SO) potential Vso-— The Schrodinger 
equation, 



;V 2 + Vp S (r) 



(1) 



is solved by using a linear combination of bulk bands 
(LCBB) method, 41 in a basis {</>„ -£ ^(k)} of Bloch or- 
bitals of band index n and wave vector k of material A 
(= InAs, InP), strained uniformly to strain e. For each 
material A, we used a basis set of n= 6 bands (including 
spin) for the hole states and n— 2 for the electron states, 
on a 6 x 6 x 16 k-point mesh around the T point. It has 
been shown that the energy levels changes in InAs/GaAs 
QDs due to the piezoelectric effects are quite small^ Be- 
cause the lattice mismatch in the InAs/InP QDs is only 
half of that of the InAs/GaAs QDs, we expect that the 
piezo-effect should be even smaller in the InAs/InP dots, 
and therefore the piezoelectric potential is neglected in 
the calculations. 

Due to the spatial confinement, the carriers in the QDs 
have strong Coulomb interactions. The many-particle 
Hamiltonian reads, 

H = J2 ^ a ,Aa,a ^M^j^U (2) 

where e is the single particle levels obtained in Eq. ([I]). 
TV, is the Coulomb integral matrix, i.e., 



- k,l ~ / j 
a, a' 



dvdr 



,^(r,a)^(r , ,(7')^(r,a)Vr(r',0 



e(r 



(3) 

where e(r — r') is the full screened dielectric constant. 57 
The many-particle Hamiltonians are solved by using 
a configuration interaction (CI) method in which the 
Hamiltonians are diagonalized in the Slater determinant 
basis^ This method has been successfully applied to 
study the optical properties of InAs/GaAs QDs, and 




FIG. 1: A sketch of the geometry of a lens-shaped InAs QD 
embedded in the InP matrix on the top of two-monolayer 
wetting layer. 
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FIG. 2: The single-particle energy levels (in meV) in a lens- 
shaped (a) InAs/InP QD and (b) InAs/GaAs QD, with di- 
ameter D= 25 nm and height h= 2.5 nm. 



the obtained results are in excellent agreement with 
experiments^ 5 -^'^^'^ 



III. RESULTS 

Figure [5] (a) and (b) show the single-particle energy 
levels for typical InAs/InP and InAs/GaAs QDs, re- 
spectively, with height h—2.5 nm and diameter D=25 
nm. The differences between the single-particle electronic 
structures of the two dots have been addressed in Ref. 
48]. Here we merely outline the main differences be- 
tween the two dots that are related to this work, (i) The 
confinement of electron and hole in InAs/GaAs QDs is 
490 and 294 meV, while in InAs/InP QDs is 325 and 
535 meV, respectively. Therefore, holes (electrons) are 



3 





9^ n - 




20.0 - 




- 

15.0 - 


r s 


10.0 - 


-mi 






5.0 - 


H — ' 




CD 


0.0 - 



(a) 



~l 1 1 1 1 

XX 



l — i — i — i — | — i — r 

InAs/lnP 



_i i i i_ 



x 

M 



_l I I L 



_L 



CD 

H — » 





50 


- 


o 






H — ' 


40 


- 


to 












cti 


30 


- 




f— 


20 


- 




10 


- 







- 



0.865 0.870 0.875 

i 1 1 1 1 1 1 1 1 1 1 1 r~ 



'(b) 



_l I I l_ 




InAs/GaAs 



_l I I l_ 



J_ 



1.025 1.030 1.035 

Transition Energy (eV) 

FIG. 3: (Color online) The PL spectra of X, XX, X+ and 
X" in a lens-shaped (a) InAs/lnP QD and (b) InAs/GaAs 
QD, with D— 25 nm and h — 2.5 nm 



much more (less) confined in InAs/lnP QDs. (ii) Due 
to the much smaller confinement potential for electrons 
in the InAs/lnP QDs, fewer electron states are confined, 
(iii) The p-orbit energy splitting of electrons and holes 
in InAs/GaAs QDs are about 1.8 meV and 7.5 meV re- 
spectively, whereas in InAs/lnP QDs, they are 1.4 meV 
and 0.7 meV, respectively, (iv) The hole states in tall 
InAs/GaAs QDs can be localized in the interface, but in 
InAs/lnP QDs, hole levels have well defined shell struc- 
ture and no localization has been observed. These dif- 
ferences can greatly influence the optical properties of 
InAs/lnP QDs, as will be shown in the following sections. 



A. Binding energies of X + , X and XX 

A typical PL spectra for lens-shaped InAs/lnP and 
InAs/GaAs QDs with D = 25 nm and h = 2.5 nm are 
presented in Fig. [3] (a) and (b), respectively. We focus on 
the alignment of emission line of \ q (x q = X, XX, X + 
and X~). The calculated alignment of x q m InAs/GaAs 
QDs agrees well with experimental results in Refs. [6ll61r- 
[fjij . and previous calculations by Narvaez et a? i 65 ' 66 . 

The binding energy of x q is defined as, 



AE b ( x i) = u( x q )-uj(X), 



(4) 



where Lo(x q ) is the corresponding recombination energy. 
In the Hartree-Fock approximation, uj{x q ) can be cal- 
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FIG. 4: (Color online) Diect Coulomb interactions Jss h \ 
and jif e ' as functions of dot height for lens-shaped (a) 
InAs/lnP QDs and (b) InAs/GaAs QDs. The diameter of the 
QDs is fixed to 25 nm. 



culated as 



65-67 
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(hh) 
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2J t 

(ee) 



(eh) 



u{X) 



j£ H \ L0(X + ) 



lo(X-) = E g + Jt ] - 2J s ( f \ oj(XX) 



J: 



(hh) 



3J ; 
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where E g = E e s — E^ (see 
Fig. [2]), j'tf' (jis h ^) is the direct Coulomb interaction 
between two electrons (holes), and Jsf^ is the direct 
Coulomb interaction between electron and hole of s lev- 
els. Here, we neglect the exchange interactions, which 
is much smaller than the direct Coulomb interactions. 
With these analytical expressions, we have 

AE b (X+)W = JW- 4f\ 

AE b (X-)^ = Jt ] - (5) 

AEtiXX)^ = j£V + - 2j(f >. 

The values of , J s ( f ) and for InAs/lnP QDs 

and InAs/GaAs QDs are shown in Fig. |4j For InAs/lnP 

QDs, we see jis^ > Jst^ > JsT^ ■ This is because 
in InAs/lnP QDs, holes are much more localized than 
electrons. In contrast, in InAs/GaAs QDs, we find the 

same order for Jss' l \ Jss h \ JsT\ for /i < 3 nm, but 
reversed order when h > 3.0 nm. Equation ([5]) can be 
used to explain the alignment of \ q of InAs/lnP QDs as 
shown in Fig. [3] (a) , but can not explain the alignment 
of x q m InAs/GaAs QDs. This is because in InAs/lnP, 

J ( S s h) - J { s e s h) - 6.2 meV, and J { s f } - J s ( s ee) - 3.2 meV, 
much larger than the correlation energies, which is about 

1 meV. However, in InAs/GaAs QDs, JsT\ jis^ and 
Jss^ are very close. So the correlation energies play very 
important roles to determine the final alignment of x q ^ 
We therefore expect that the alignments (and binding 
energies) of Xq ln InAs/lnP QDs are very different from 
those of InAs/GaAs QDs. 

We show the binding energies of X + , X~ and XX 
as functions of dot height, dot diameters and aspect ra- 
tio (ano/ai-io) in Fig. [5] (a), (b) (c) respectively. We 
find that the binding energies of X + and XX decrease 
dramatically with the increasing of dot height, whereas 
the binding energy of X~ increases. For a lens-shaped 
dot with diameter D=20 nm, the binding energy of XX 
equal zero, i.e., lu(XX) = oj{X) at dot height h ~ 4nm. 
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FIG. 5: (Color online) (a) The binding energies as functions 
of dot height for lens-shaped QDs with D—20 nm. (b) The 
binding energies as functions of dot diameter for lens-shaped 
QDs with h — 3 nm, (c) The binding energies for elongated 
QDs, where the total volume is fixed to 400 nm 2 , and h=3.5 
nm. 



At this special point, the dot can be used as an entangled 
photon emitter via the time-reordering scheme proposed 
by Avron et al£& In contrast, the binding energy of X + is 
always positive, whereas the binding energy of X~ is al- 
ways negative, even though they all become small for very 
tall dots. As shown in Fig. 0(c), there is no obvious shift 
of binding energies in elongated QDs when dno/ai-io 
changes from 0.9 to 1.1. The observed trends of the bind- 
ing energies for InAs/InP QDs differ greatly from those 
of InAs/GaAs QDs. For example, in InAs/GaAs QDs, 
the binding energies of X + and X~ can be either posi- 
tive or negative^ due to the cross over of the Coulomb 
interactions shown in Fig. @] . These differences can be 
traced back to the difference of confinement potentials 
between the two dots. 



B. Hidden correlation 
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FIG. 6: (Color online) The hidden correlations as functions 
of (a) dot height, (b) dot diameter, and (c) aspect ratio in 
InAs/InP QDs. 



consistent Hartree-Fock approximation, the binding en- 
ergies of XX, X + and X~ cancel each other exactly 
(see also Eq. fj5j), i.e., A HF = 0. Therefore the value 
that A ^ arise purely from the correlation effect. In 
Ref. 69], it is shown that A varies between 1.2 - 2.2 
meV, for a wide range of geometry and exciton energies 
of InAs/GaAs QDs, which agrees well with all available 
experimental data from different groups^ The hidden 
correlation is also confirmed by the effective mass quan- 
tum Monte Carlo simulations.— 

The hidden correlations in InAs/InP QDs is shown in 
Fig. HI as functions of dot heights, diameters and as- 
pect rations. We find A = 0.7 ~ 0.9 meV for all QDs 
with different sizes and geometries studied here. The hid- 
den correlation of InAs/InP QDs is smaller than that of 
InAs/GaAs QDs, which can be understood as follows. In 
InAs/InP QDs, the hole is more confined, so the energy 
level spacings are much larger than those in InAs/GaAs 
QDs (see Fig. [2]) and therefore the correlation energies 
are reduced (see Ref. Experimental confirmation 

is called for this prediction. 



The exciton binding energies change dramatically with 
the geometry and compositions of QDs. Recently pro- 
posed by some authors of the paper that the hidden cor- 
relation A of QDs are always positive and remarkably 
constant for a large range of self-assembled QDs as a 
consequence of the Coulomb correlation effects^ The 
hidden correlation is defined as, 

A = AE b (XX)-AE b (X-)-AE b (X+), (6) 

where AE b (x q ) nas been defined in Eq. Q and can 
be extracted directly from Fig. [S] Under the non-self- 



C. Lifetime of exciton, biexciton and trions 

The radiative recombination lifetimes of exciton, biex- 
citon and trion are very important parameters for many 
applications of the QDs. At low temperature, the lifetime 
is mainly determined by the dipole matrix element, 

M id = (f\n-p\i), (7) 

where p is the dipole momentum and n is the polarization 
vector of the electromagnetic field. The recombination 
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FIG. 7: The radiative lifetimes of (a) X and (b) XX in the 
lens-shaped InAs/InP QDs as functions of dot height. 
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FIG. 8: The radiative lifetimes of (a) X~ and (b) X + in the 
lens-shaped InAs/InP QDs as functions of dot height. 



lifetime r t t at zero temperature can be calculated from, 

- = fm\M,A\ (8) 

where n r is the reflective index, and — \u>i — u>t\ is 
the recombination energy, mo is the mass of the elec- 
tron and c is the velocity of light in vacuum. The lin- 
ear dependence of 1 /tj / on refractive index is applicable 
only when the QD and the matrix has similar dielectric 
constantS i 66 i 70 At finite temperature, the lifetime can be 
calculated as^ 

1 v - ^ ni 

where rij is occupation number of the initio state accord- 
ing to the Boltzmann distribution, and m = 1. As- 
suming that the thermalization between the dark and 
bright states is much longer than the exciton life time, 
we only take the bright states into consideration, which 
is somehow different from the treatment in Ref. [661 ] . 

The radiative lifetimes of exciton in lens-shaped 
InAs/InP and InAs/GaAs QDs are compared in Fig. [Jj 
(a). The major differences between the two dots are: (i) 
The exciton lifetime in InAs/InP is almost twice longer 
than that of InAs/GaAs for flat QDs. There are two 



reasons for the longer lifetime in InAs/InP QDs. First, 
InAs/InP QDs have smaller exciton energies and second, 
InAs/InP QDs have much smaller dipole moment, be- 
cause the electrons are less confined. Experimentally, the 
lifetimes of X and XX in InAs/GaAs QDs are found to 
be around 0.8 - 1.2 ns^^ ' 16 ' 18 i 63 very close to the calcu- 
lated results presented in Fig. [JJ The measured exciton 
lifetimes in InAs/InP QDs by Sasakura et alM^ are about 
2.0 ns, which also agree very well with the calculated val- 
ues in this work, (ii) The exciton lifetime in InAs/InP 
QDs decreases with increasing of the dot height. In the 
InAs/GaAs QDs, the opposite trend is found, because the 
hole wave functions localize more on the interface of the 
dots4£ (hi) For flat QDs (h < 3.5 nm), we find that the 
exciton lifetime in InAs / GaAs dots is almost independent 
of the dot diameter. However, in InAs/InP QDs, the ex- 
citon lifetime depends strongly on the dot diameter. For 
example, when D = 20 nm, h= 2.5 nm, t(X)= 2.26 ns. 
When we increase D to 25 nm, t(X) decrease to 2.04 ns. 
This is also due to the weak confinement potential for 
electrons in the InAs/InP QDs. Increasing the diameter 
of a QD results in more confined electron state, and thus 
lager overlap between electron and hole wave functions 
and lager dipole moment. 

We further compare the lifetimes of biexciton in Fig. [7] 
(b), and of charged exciton between the two dots in Fig. 
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[8j The values of lifetime for these exciton complexes 
are very different. However the trends of lifetimes as 
functions of dot geometry are rather similar. Roughly 
speaking, the charged excitons (X + , X + ) and X have 
lifetimes about twice longer than the biexciton^ i.e., 

t{X+) ~ t(X~) ~ t(X) ~ 2t(XX). (10) 



D. Phase diagrams of electrons and holes 

The studies of electron/hole filling process in self- 
assembled QDs are stimulated by the experiments that 
allow the electron/hole be loaded into the QDs one by 
one and measure the charging energiesi 19 i 71 Unlike the 
atoms, where large Coulomb interaction J ~ 10 eV en- 
sure the Hund's rule (maximum spin) and Aufbau prin- 
ciple. In self-assembled QDs, the Coulomb interaction is 
about 20 - 25 meV for electron-electron interaction and 
15 - 20 meV for hole-hole interaction, which is the same 
order of the energy level spacings in QDs. So, the elec- 
tron/hole filling process should be quite different from 
that in atoms. In strongly confined QDs (dimeter ~ 3 
nm), it is found that Hund's rule is generally obeyed, 
but violation of the Aufbau principle is a common fea- 
ture when single-particle energy levels separated by a 
few meVi££ In the self-assembled InAs/GaAs QDs, He 
et a f 2 ^ 73 found that it is possible for holes to violate 
both the Hund's Rule and the Aufbau Principe, which is 
confirmed by the PL spectra of highly charged excitons 
X N + and X N - by Ediger et al^.. 

Let us first look at the charging energy and addi- 
tion energy of the QDs, which can be measured di- 
rectly in the experiments. The charging energy n(N) = 
E(N) — E(N — 1), measures the energy needed to load 
an additional carrier into the QDs. The addition energy 
is the energy difference between fJ,(N) and fi(N — 1), i.e., 
A (AT — 1, N) = n(N) - fi(N - 1) = E(N) - 2E(N - 1) + 
E(N — 2). In the Hartree-Fock approximation, the total 
energy of N electrons confined in QDs can be calculated 



as 



57 



E(N) 



i<j 



M^i + > X J i,3 



)nim 



(11) 



where rij is the occupation number of the z-th level. The 
diagonal Coulomb energy J,- ■ — r ' ' 



'h3 



Vl'l and the exchange 



energy K hj = Lyj (see Eq. [3]). 

The electron/hole addition energies are presented in 
Fig. l(a),(c), for InAs/GaAs dots and in Fig. 1(b) , (d) 
for InAs/InP QDs, respectively. The electron addition 
energies of InAs/InP QDs are much smaller than those 
of InAs/GaAs QDs of similar sizes. However, the hole 
addition energies of InAs/InP QDs are much larger than 
those of InAs/GaAs QDs due to their much larger single- 
particle level spacings, as illustrated in Fig. [2] More in- 
terestingly, where the hole addition energy in InAs/GaAs 
QDs show two peaks at N=3 and 5, the hole addition en- 
ergy in InAs/InP QDs show only a strong peak at N—3, 



which is similar to that of electrons. The reason for this 
difference is that the hole energy levels have well defined 
shell structure in InAs/InP QDs, but not in InAs/GaAs 
QDs. 

We further look at the filling orders of the electron/hole 
in the InAs/InP QDs. The filling orders of electrons and 
holes in QDs are determined by the following factorsi 72 i 73 
(i) the single-particle energy level spacing, (ii) the p, 
d energy level splittings and (iii) the Coulomb interac- 
tions. We adopt the general approaches in Refs. [zHzl], 
and use the reduced single particle energy level spac- 
ing 5{p 2 ,Pi)/J ss = \E P2 - E P1 \/J SS and 5{d 1 ,p 2 ) / J ss = 
^ E P2 \/J SS to construct the phase diagrams of few 
carriers in QDs. Experimentally, the p-orbit splitting 
8{p2 1 pi) in QDs can be measured via capacitance- voltage 
spectroscopy^ or terahertz pump-probe method^ 

For N — 4, 5, 6, the possible ground state configu- 
rations represented by the spectroscopic notation read 



as 



73 
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N = 5 



N = 6 
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4 A: 
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3 n = 
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i E * = 



«V)(pI)(p!), 

} V)(rfpt); 

'V)(rfx4)(4), 
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.t«4. 



)(p1)(pI)(4)(4 



(12) 



(13) 



( s V)(pMx44). 



(14) 



phase diagrams of electrons 



The phase diagrams of electrons in InAs/InP QDs 
are presented in Fig. [TOTb). contrast to those from 
two dimensional effective mass approximation (2D-EMA) 
calculations. The phase diagrams are determined by 
S(p2,Pi)/ J S s and 5(d\,pz)/ J aa . Previous calculations 48 
have shown that the electron levels in InAs/InP QDs have 
well defined electron structure, similar to those given by 
2D-EMA. Therefore the phase diagram from 2D-EMA 
and atomistic model share quite similar features. The 
ground state configurations for N — 4 - 6 are 3 S, 2 n, 4 E 
(indicated by the red spots in the figure), respectively, fol- 
lowing Hund's rule and Aufbau principle. However, for 
A^=4, the area allow stable 3 E configuration in InAs/InP 
QDs is much smaller than that of the 2D-EMA model. 



2. phase diagrams of holes 

The phase diagrams of holes in InAs/InP QDs are pre- 
sented in Fig. [TOTc). Unlike those of the InAs/GaAs 
QDs, the hole levels in InAs/InP QDs have well defined 
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FIG. 9: (Color online) The electron addition energies for (a) 
InAs/GaAs QDs and (b) InAs/lnP QDs. The hole addition 
energies are shown in (c) for InAs/GaAs QDs and in (d) for 
InAs/lnP QDs. 



shell structure. The ground state of N = 3 is 3 E, satis- 
fying Hund's rule. However, the phase diagrams of holes 
in the InAs/lnP QDs still differ greatly than calculated 
by 2D-EMA for N = 5, 6. 

N = 5: The are three possible ground state configura- 
tions for 2D-EMA model and in InAs/GaAs QDsiS He et 
have predicated that the non-Aufbau phase 2 A can 
be found for large 5(p2,Pi) in InAs/GaAs QDs. Similar 
phase, however, can not be found in InAs/lnP QDs. This 
can be understood as following. In InAs/lnP QDs, the p- 
orbital splitting 5(p2,Pi) ~ 1-4 meV and p-d level spac- 
ing 8{d\,p2) ~ 15 - 30 meV (see Fig. [21), so the Coulomb 
interactions between holes are not strong enough to over- 
come the p-d spacing in InAs/lnP QDs. Therefore 2 A is 
always unfavorable in energy in InAs/lnP dots. 

N = 6: There are four possible ground state configura- 
tions for TV = 6. In both InAs/GaAs and InAs/lnP QDs, 
the 5 A can not be observed. However, there are still ob- 
vious difference between the filling orders in InAs/GaAs 
and InAs/lnP QDs. In InAs/GaAs QDs, 75 the hole 
can be localized in the interface, resulting in extremely 
small level spacing between pi and d\ orbits, hence the 
which violate the Hund's rule, is possible! 52 ' 72 ^ 3 
In InAs/lnP QDs, only two phases, : S and 3 n, can be 
found. 

The hole energy level spacings in InAs/lnP QDs are 
significantly larger than in InAs/GaAs dots, and the p 
level splittings are significantly smaller. Therefore, the 
ground states for N=5, 6 in InAs/lnP QDs are 2 n and 
1 S respectively, as indicated by the red dots in FigEJc), 
which obey both Hund's rule and Aufbau principle. 
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FIG. 10: Phase diagrams of InAs/lnP QDs with number of 
electrons/holes N = 4, 5, 6. (a) The results from 2D-EMA 
calculations for both electrons and holes, which are identi- 
cal for the parabolic confinement potential, (b) The results 
of InAs/lnP QDs from atomistic pseudopotential method for 
electrons and (c) for holes. The red dots indicate the config- 
urations of the ground state. 



E. PL Spectra of highly charged excitons 

Having clarified the ground state configurations and 
filling orders of electrons and holes in the InAs/lnP QDs, 
we now turn to the PL spectra of highly charged excitons, 
which provides a useful tool to explore the complex inter- 
actions in the QDs.— In InAs/GaAs QDs, it is found that 
the PL spectra of highly charged excitons exhibit some 
peculiar properties due to the breakdown of the Aufbau 
principle for hole^ 2 - 

We present the calculated PL spectra of X N+ and 
X N - [N = 0-5) for InAs/lnP QDs in Fig. [TTJ We also 
calculated the PL spectra of highly charged InAs/GaAs 
QDs. The results (not shown) agree very well with pre- 
vious calculations in Ref. [52[. There are several dif- 
ferences between the PL spectra of these InAs/lnP and 
InAs/GaAs dots. 

(i) The transition intensities of X N+ and X N ~ in 
InAs/lnP are weaker than that in InAs/GaAs QDs, 
due to the relatively longer lifetime in InAs/lnP QDs. 
Rough estimation shows that the transition intensities in 
InAs/lnP are half of those in InAs/GaAs QDs. 

(ii) In InAs/lnP QDs, the main peak of X 5 ~ is on the 
lower energy side of X° by about 16 meV, and that of 
X 5+ is on the higher energy side of X° by about 12 meV, 
much larger than the energy shifts in InAs/GaAs QDs, 
where both X 5 ~ and X 5+ are on the lower energy side 
of X° by about 10 meV and 2 meV respectively. These 
differences are due to the different confinement potentials 
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(iv) In the InAs/GaAs QDs, holes may violate Hund's 
rule and Afubau's principle. For example, in the 
InAs/GaAs QDs, X 3+ has a close shell ground state,— 
and therefore no polarization dependent photolumines- 
cence. In contrast, X 3+ of InAs/InP QDs has a open 
shell ground state, and therefore has polarization depen- 
dent photoluminescence. Furthermore, there are many 
more transition peaks due to the spin splitting in the 
InAs/InP spectra than in the InAs/GaAs spectra. Pre- 
vious calculations^ show that X 5+ has open shell ground 
state and has polarization dependent photoluminescence 
in the InAs/GaAs QDs. However, in the InAs/InP QDs, 
X 5+ has close shell ground state, and no polarization 
dependence in the photoluminescence. These differences 
reflect that the hole occupation in the two dots are sig- 
nificantly different. 



IV. SUMMARY 



FIG. 11: The PL spectra of highly charged excitons in a lens- 
shaped InAs/InP QD with D=25 nm and h = 2.5 nm at 4 
K. The number in each panel is the pecentage of the leading 
configuration of the initial state CI wave functions. 



of the two dots, as discussed in Sec. IIII Al 

(iii) The number in each panel is the overlap between 
the many-particle wave functions from CI calculations 
and those from Hartree-Fock approximations for the ini- 
tial states, i.e., 



0(X a ) = \(* cl (X a )\V HF (X a )}f 



(15) 



0(X a ) represents the percentage of the leading config- 
urations in state X a . We find that there are signifi- 
cant differences of 0(X a ) between X N ~ and X N+ in 
InAs/InP QDs. For X N ~, O is generally larger than 
98%, in contrast to 95% in InAs/GaAs QDs. More- 
over, the configuration hybridization of the final states 
Q £ x N ~ (N — 4, 5) is also much smaller than that in 
the InAs/GaAs QDs. Therefore, A 4- and X 5 " have 
much fewer peaks in the InAs/InP PL spectra than in 
the InAs/GaAs spectral. For X N+ , we find that the 
hybridization effect in InAs/InP is much stronger than 
that in InAs/GaAs QDs. For example, in InAs/InP QDs, 
0(X 2+ ) = 56%, much smaller than 84% in InAs/GaAs 
QDs. This giant difference is from the smaller p level en- 
ergy splitting in InAs/InP QDs, as illustrated in Fig. [2] 
The strong hybridization effects lead to the much more 
complex X N+ PL spectra in InAs/InP QDs when N > 
3. 



In this work, we present a comprehensive studies of 
the optical properties of InAs/InP QDs using an empiri- 
cal pseudopotential method and configuration interaction 
treatment of many-particle interactions. The results are 
contrast to those of InAs/GaAs QDs, whose properties 
are well understood, both theoretically and experimen- 
tally. The main difference between the optical properties 
of the two QDs can be summarized as: (i) The align- 
ment of emission lines of neutral exciton, charged exci- 
ton and biexciton in InAs/InP QDs is quite different from 
that in InAs/GaAs QDs. (ii) The hidden correlation in 
InAs/InP QDs is about 0.7 - 0.9 meV, much smaller than 
that in InAs/GaAs QDs (1.2 - 2.2 meV). (iii) The life- 
times of neutral exciton, charged exciton and biexciton 
in InAs/InP QDs are about twice longer than those in 
InAs/GaAs QDs. (iv) The phase diagrams of electrons 
and holes in InAs/InP QDs differ greatly from those in 
InAs/GaAs QDs. The filling orders of electrons and holes 
are shown to obey the Hund's rule and Aufbau principle, 
(vi) The PL spectra of highly charged excitons show some 
significant difference between the two dots, due to the 
different filling orders and Coulomb interactions between 
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